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Abstract. In recent years two nonlinear dispersive partial differential equa- 
tions have attracted a lot of attention due to their integrable structure. We 
prove that both equations arise in the modeling of the propagation of shal- 
low water waves over a flat bed. The equations capture stronger nonlinear 
effects than the classical nonlinear dispersive Benjamin-Bona-Mahoney and 
Korteweg-de Vries equations. In particular, they accomodate wave breaking 
phenomena. 



1. Introduction 

The study of water waves is a fascinating subject because the phenomena are 
familiar and the mathematical problems are various cf. [33] . Due to the relative 
intractability of the governing equations for water waves in regard to inferring from 
their direct study qualitative or quantitative conclusions about the propagation of 
waves at the water's surface, from the earliest days in the development of hydrody- 
namics many competing models were suggested. Until the second half of the 20'th 
century, the study of water waves was confined almost exclusively to linear theory 
[T2] . While linearisation gives insight for small perturbations on water initially at 
rest, its applicability fails for waves that are not small perturbations of a flat water 
surface. For example, linear water wave theory gives no insight into the study of 
phenomena which are manifestations of genuine nonlinear behaviour, like breaking 
waves breaking and solitary waves [31j . Many nonlinear models for water waves 
have been suggested to capture the existence of solitary water waves and the asso- 
ciated phenomenon of soliton manifestation |21j . The most prominent example is 
the Korteweg-de Vries (KdV) equation [23], the only member of the wider family 
of BBM-type equations [5] that is integrable and relevant for the phenomenon of 
soliton manifestation [15]. Another development of models for water waves was 
initiated in order to gain insight into wave breaking, one of the most fundamental 
aspects of water waves for which there appears to be no satisfactory mathematical 
theory [33]. Starting from the observation that the strong dispersive effect incor- 
porated into the KdV model prevents wave breaking, Whitham (see the discussion 
in [33| ) initiated the quest for equations that are simpler than the governing equa- 
tions for water waves and which could model breaking waves. The physical validity 
of the first proposed models is questionable but two recently derived nonlinear 
integrable equations, the Camassa-Holm equation [7] and the Degaspcris-Procesi 
equation [14] , possess smooth solutions that develop singularities in finite time via a 
process that captures the essential features of breaking waves cf. [33] : the solution 
remains bounded but its slope becomes unbounded. Our aim is to prove the rele- 
vance of these two equations as models for the propagation of shallow water waves, 
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proving that both are valid approximations to the governing equations for water 
waves. In our investigation we put earlier (formal) asymptotic procedures due to 
Johnson |22j on a firm and mathematically rigorous basis. We also investigate in 
what sense these two models give us insight into the wave breaking phenomenon 
by some simple numerical computations. 

1.1. Unidirectional asymptotics for water waves. For one dimensional sur- 
faces, the water waves equations read, in nondimensionalized form, 

![xd%3> + d z $ 2 = in O t , 

d z § = 0, at z = 0, 

d t C-j i (-nd a! (d x $ + d 1 .*) = at z = e(, 

fit* + i(0**) a + ^(0,$) a = at z = e(, 

where x i— > s£(t, x) parameterizes the elevation of the free surface at time t, fit = 
{(x, z), —1 < z < e((t, x)} is the fluid domain delimited by the free surface and the 
flat bottom {z = —1}, and where $>(t, ■) (defined on £l t ) is the velocity potential 
associated to the flow (that is, the two-dimensional velocity field v is given by 
v = (9 X $, d z $) T ). Finally, e and [i are two dimcnsionless parameters defined as 

a h 2 
£= h> »=V> 

where h is the mean depth, a is the typical amplitude and A the typical wavelength 
of the waves under consideration. Making assumptions on the respective size of e 
and fj,, one is led to derive (simpler) asymptotic models from ([I]). 
In the shallow- water scaling (/i <C 1), one can derive the so-called Green-Naghdi 
equations (see [19] for the derivation, and [3j for a rigorous justification), without 
any assumption on e (that is, e = 0(1)). For one dimensional surfaces and flat bot- 
toms, these equations couple the free surface elevation £ to the vertically averaged 
horizontal component of the velocity, 
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(2) u(t, x) = - r - I7 I d x <f>(t,x, z)dz; 



(3) 



1 + eC 
and can be written as 

Ct+[(l + eC)«L = 

u t + ( x + suu x = | xq^[(l + £() 3 (u xt + euu xx - eu 2 x )] x , 
where 0(p 2 ) terms have been discarded. 

If we make the additionnal assumption that £<1, then the above system reduces at 
first order to a wave equation of speed ±1 and any perturbation of the surface splits 
up into two components moving in opposite directions. A natural issue is therefore 
to describe more accurately the motion of these two "unidirectional" waves. In the 
so-called long- wave regime 

(4) /i«l, e = 0(n), 

Korteweg and de Vries [55] found that, say, the right-going wave should satisfy the 
KdV equation 

3 1 
ut +u x + e-uu x + (J,-u xxx = 
2 b 

(and C = u + 0(s,fi)), which at leading order reduces to the expected transport 
equation at speed 1. More recently, it has been noticed by Benjamin, Bona, Ma- 
honey [5] that the KdV equation belongs to a wider class of equations (the BBM 
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equations, first used by Peregrine [29j and sometimes also called the regularized 
long-wave equations) which provide an approximation of the exact water waves 
equations of the same accuracy as the KdV equation: 

3 1 

(5) u t + u x + -euu x + fi(au xxx + (3u xxt ) = 0, with a -(i=-. 

I o 

The equations ((5|) contain both nonlinear effects (the uu x term) and dispersive 
effects (the u xxx a,mi u xx t terms) due to the scaling (|4|). However, these equations do 
not account correctly for large amplitude waves, whose behavior is more nonlinear 
than dispersive. For such waves, characterized by larger values of e, it is natural to 
investigate the following scaling (which we call Camassa-Holm scaling): 

(6) n « 1, e = 0(Vm). 

With this scaling, one still has e < 1 and thus the same reduction to a simple 
wave equation at leading order; the dimensionlcss parameter is however larger here 
than in the long wave scaling, and the nonlinear effects are therefore stronger. In 
particular, a stronger nonlinearity could allow the appearance of breaking waves - a 
fundamental phenomenon in the theory of water waves that is not captured by the 
BBM equations. We show in this paper that the correct generalization of the BBM 
equations ((5]) under the scaling (|6|) is provided by the following class of equations: 
3 

(7) u t + u x + -euu x + ix{au xxx + 0u xxt ) = s^uil xxx + Su x u xx ) 
(with some conditions on a, /?, 7 and S). 

Notice that for an equation of the family ([7]) to be well-posed it is necessary that 
(3 < 0, as one can see by analyzing the linear part via Fourier transforms. 

We want to insist on the fact that ([7]) provides an approximation of the same 
order 0(/x 2 ) as the BBM equations (J5|) to the Green-Naghdi equations. The only 
difference in the derivation of these equations lies in the different scalings (j4|) and 
([6|). Since of course O(^l) = 0( v / 7i) when /j, is small, the long-wave scaling ((4]) 
is contained in the CH scaling ©, and consequently, the BBM equations can be 
recovered as a specialization of ([Jj when e = O(fi) and not only 0(y/JI). 

1.2. The Camassa-Holm and Degasperis-Procesi equations. Among the 
various type of equations ([7]) with < there are only two with a bi-Hamiltonian 
structure: the Camassa-Holm and the Degasperis-Procesi equations [2D]. Notice 
that while the KdV equation has a bi-Hamiltonian structure (see |15j). this is not 
the case for the other members of the BBM family of equations ([5]). The impor- 
tance of a bi-Hamiltonian structure lies in the fact that in general it represents 
the hallmark of a completely intcgrable Hamiltonian system whose solitary wave 
solutions are solitons, that is, localized waves that recover their shape and speed 
after interacting nonlincarly with another wave of the same type (see |15[ I21| ). 

1.2.1. Camassa-Holm equations. The Camassa-Holm (CH) equations are usually 
written under the form 

(8) U t + kU x + WU X - U txx = 2U X U XX + UU XXX , 

with k € M. A straightforward scaling argument shows that if k ^ 0, ||SJ can 
be written under the form (J7j by setting u(t,x) = aU{b{x — vt),ct) and a = 
b 2 = —^7, v = jj, c = i(l — v) (which requires (i < and leads to 7 — — ^ and 
S = 27). This motivates the following definition: 
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Definition 1. We say that £?y is a C amass a- Holm equation if the following con- 
ditions hold: 

/3<0, a?0, p = -2ry, S = 2ry. 
For all k 7^ 0, the solution u to ([?]) is transformed into a solution U to (0) by the 
transformation 

1 ,x v t, 
U(t,x) = -u(- + -t,-), 
ao c c 

with a = -4(1 — v), b 2 = —4-, v = §, and c = k(l — v). 

First derived as a bi-Hamiltonian system by Fokas & Fuchssteiner [17j , the equa- 
tion ((SJ) gained prominence after Camassa-Holm [7] independently re-derived it as 
an approximation to the Euler equations of hydrodynamics and discovered a num- 
ber of the intriguing properties of this equation. In [7] a Lax pair formulation of 
© was found, a fact which lies at the core of showing via direct and inverse scat- 
tering [SI [TU] that ([5]) is a completely integrable Hamiltonian system: for a large 
class of initial data, solving (JSJ) amounts to integrating an infinite number of lin- 
ear first-order ordinary differential equations which describe the evolution in time 
of the action-angle variables. The Camassa-Holm equations shares with KdV this 
integrability property as well as the fact that its solitary waves are solitons [TTJ [TUJ . 
We refer to [TH] for a discussion of these properties in the context of the KdV model. 

1.2.2. Degasperis-Procesi equations. The Degasperis-Procesi (DP) equations are 
usually written under the form 

(9) U t + W x + 4UU X - U txx = W X U XX + UU XXX , 

with k G M. The same scaling arguments as for the CH equation motivate the 
following definition: 

Definition 2. We say that fiTty is a Degasperis-Procesi equation if the following 
conditions hold: 

/3<0, /3 = -^7 5= 3 7 . 

For all k ^ 0, the solution u to ^ is transformed into a solution U to (0) by the 
transformation 

U(t,x) = -u(- + -t,-), 
a b c c 

with a = ^(l - v), b 2 = v=% and c = |(1 - v). 

Equation ([9|) , first derived in |14j , is also known to have a Lax pair formulation 
[T5] and its solitary waves interact like solitons [55] . Just like the KdV equation (see 
[15] ) and the Camassa-Holm equation (see [24]), the Degasperis-Procesi equation 
has infinitely many integrals of motion. 

1.3. Wave breaking. In addition to the properties of ([8]) and ([9j mentioned be- 
fore, the importance of these two equations is enhanced by their relevance to the 
modeling of wave breaking, one of the most important but mathematically still 
quite elusive phenomena encountered in the study of water waves. 

Definition 3. We say that there is wave breaking for an equation of the form Q), 
if there exists a time < t 6 '^ < oo and solutions u to ^ such that 

u G i°°([0, t s,tl ] X K) and lim^ \d x u(t, -)|oo = oo. 
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The rationale of this definition is that if the flow velocity u is to first order a 
good approximation to the surface profile it is then reasonable to expect that 
the above blow-up pattern has a similar counterpart in terms of Q. If this were the 
case, the boundedncss of the wave height in combination with an unbounded slope 
captures the main features of a breaking wave |33j . 

For KdV in particular, as well as for any other member of the BBM family ((5]), 
all smooth initial data u(0, •) decaying at infinity develop into solutions defined for 
all times (see e.g. [35]) so that the BBM family [3] does not model wave breaking 
[H [30] . To remedy this shortcoming of the KdV equation Whitham proposed to 
formally replace the dispersive term u xxx by a convolution with a singular function 
chosen so that the newly obtained equation presents wave breaking (see [T8] I33j). 
However, this formal process destroys the integrability and soliton features of the 
KdV equation. In contrast to this, both ([8]) and ([9]) admit breaking waves in the 
sense of Definition [3J cf . [3 M W\ HE] for © and [S] for ©. In this paper we 
explore the wave breaking phenomenon for both equations not in the restricted 
sense provided by Definition [3j but by studying the nonlinear equation describing 
the evolution in time of the free surface. While our results vindicate the fact that 
it is appropriate to use in this context Definition [3] to describe breaking waves, 
there is a slight twist. To be more precise, let us distinguish between two types of 
breaking waves that can be observed. In a plunging breaker the slope of the wave 
approaches — oo at the breaking location as we reach breaking time, while in a 
surging breaker the slope becomes +oo. Considering the case of the Camassa-Holm 
equation ([8]), due to the fact that singularities in a smooth solution can appear only 
if infjjgR {u x (t, x)} — ► —oo as we approach breaking time while sup^gg {|u(t, x)\} re- 
mains uniformly bounded (see [5]), we would expect to observe a plunging breaker at 
the free surface. However, as we recall here, the Camassa-Holm equation describes 
the behavior of the vertically averaged horizontal component of the velocity it; the 
free surface elevation £ can be given in terms of u = u + ju 2 + ^u xt + O(e^t), 
see Proposition [1] below) but such an asymptotic expression of course breaks down 
when u x becomes singular, and cannot be used to describe the behavior of the free 
surface when there is "wave breaking" for the velocity. Since "wave breaking" is 
a very intuitive notion when it refers to the free surface elevation, we show in this 
article that it is possible to revert the usual approach; that is, we derive an evo- 
lution equation for the surface elevation £ and give an asymptotic expression for 
the velocity u in terms of C (cf . £|2.2p . This allows us to prove that wave breaking 
indeed occurs for the surface elevation but that, as opposed to what happens for 
the velocity, this is a surging breaker! The difference between plunging and surging 
breakers is graphically illustrated by numerical computations in N3.4I 



2. Derivation of asymptotical equations for the unidirectional limit 



We derive here asymptotical equations to the Green-Naghdi equations in the 
Camassa-Holm scaling ([6]). We recall that the Green-Naghdi equations are given 
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Since we work under the Camassa-Holm scaling, we restrict our attention to values 
of e and /i satisfying 

(10) (e,n) EV:= {n E (0,Ho),e< M^fji}, 

for some p,Q > and M > 0. Equations for the velocity u (including the CH 
and DP equations) are first derived in £|2.1[ and equations for the surface elevation 
£ arc obtained in i)2.2l The considerations we make on the derivation of these 
equations are related to the approach initiated by Johnson [22) . approach that 
is substantiated and extended by our analysis. In addition, we explore the wave 
breaking phenomenon. 

2.1. Equations on the velocity. At leading order, the Green- Naghdi equations 
degenerate into a simple wave equation of speeds ±1; including the 0{e) terms, one 
can easily check that the (say) right-going component of the wave must satisfy 

3 

(11) U t + U x + -EUUx = 0, 

and ( = u + 0(e). 

If we want to find an asymptotic at order 0(/j?) (recall that e = 0(fi 1 ^ 2 )), it is 
therefore natural to look for u as a solution of a perturbation of (jTTJ) including terms 
of order 0(fj,) and 0(efi) similar to those present in ([3]). Thus we want u to solve 
an equation of the form (|7|) where a, (3, 7 and 5 arc coefficients to be determined. 
We prove in this section that under certain conditions on the coefficients, one can 
associate to the solutions of a family of approximate solutions consistent with 
the Green-Naghdi equations ([3]) in the following sense: 



Definition 4. Let fi > 0, M > 0, T > and V be as defined in uuj). 

A family (C £ ''\ ^'^(e.^ev * s consistent (of order s > and on [0, ■£■]) with the 

Green-Naghdi equations if for all (e,fi) € V, 



G+ [(1 + eCH^ A!' A 



u t + Cx + £uu x = + e() 3 (u xt + euu xx - eul)] x + [i 2 r E 2 11 ] 

with (rl'i* ,r e ^) {Et ^ v bounded in L°°([0, | ], H S (R) 2 ). 

The following proposition shows that there is a one parameter family of equations 
of the form ||7J) consistent with the Green-Naghdi equations. 

Proposition 1. Let p £ K and asssume that 

1 3 1 . 9 23 

a = p. P=p—-, 7 = — — 7:, o = — —p— —. 

6 ' 2 F 6 2 F 24 

Then there exists D > such that: 

• For all s > and T > 0, 

• For all bounded family (u £ '^) (e ^ )( , v £ C([0, j];H s+D (R)) solving ffi), 

the family (£ e ' M , u e,il )(e,^)ev j with (omitting the indexes s,fi) 

el 15 
C := u + -u 2 + fi-u xt - en [-uu xx + — u 2 x ] , 

is consistent (of order s and on [0, —]) with the Green-Naghdi equations (QJ). 
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Remark 1. i. One can recover the equations (26a) and (26b) of |22] with p = — 
( and thus a = — j^ , /? = — i ,7 = — -^ and 6 = — and p = | f and thus a = 
(3 = 0, 7 = — j2 and S = — ^) respectively. 

ii. The one parameter family of equations considered in the proposition admits 
only one representant for which 8 = 27 (obtained for p = — jz, and thus 7 = 57}, 
5 = j^). However, since (3 = — X 7^ —27, the corresponding equation is not a 
Camassa-Holm equation in the sense of Definition^ 

iii. There is no possible choice of p such that 6 — 3"f in Proposition QJ Con- 
sequently, none of this one parameter family of equations is a Degasperis-Procesi 
equation. Notice that among all equations with (3 < there are only two with a 
bi-Hamiltonian structure: the Camassa-Holm and the Degasperis-Procesi equations 

W- 

Proof. For the sake of simplicity, we use the notation O(fi), 0(/i 2 ), etc., without 
explicit mention to the functional normed space to which we refer. A precise state- 
ment has been given in Definition it would be straightforward but quite heavy 
to maintain this formalism throughout the proof. 
Step 1. If it solves ([7]) then one also has 
3 

(12) u t + u x + e-uu x + [J,au xxt = efi[buu xx + cul] x + 0(/i 2 ), 

with a — (3 — a, b — 7 + |a and c = ^{5 + 3a — 7). 
Differentiating ([7]) twice with respect to x, one gets indeed 

3 

u xxx = -u xxt - -ed 2 x {uu x ) + O(fj), 

and we can replace the u xxx term of (J7J) by this expression to get (|12D . 

Step 2. We seek v such that if ( = u + ev and u solves ([7]) then the second equation 

of © is satisfied ut to a 0(/i 2 ) term. This is equivalent to checking that 

u t + [u + ev] x + euu x - -u xxt = ^- [ - uu xxt + [3uu xt + uu xx - u 2 ] x ] + 0(^ 2 ), 

= -—[uu xx + -u x ] x + 0(n) 

the last line being a consequence of the identity u t = —u x + 0(e) provided by (fT2")l . 
The above equation can be recast under the form: 

3 

£v x + [u t + u x + e-uu x + iiau xxt - e^[buu xx + cuj J 

el 11 
= -uu x + n{a + -)u xxt - £fi[(b + -)uu xx + (c + -)u x ] x + 0{^ 2 ). 

From Step 1, we know that the term between brackets in the lhs of this equation 
is of order 0(fi 2 ), so that that the second equation of ([3]) is satisfied up to 0(/i 2 ) 
terms if 



= -uu x + fi(a + -)u xxt - en[(b+ -)uu xx + (c + -)*4] x + 0{/j. 2 ), 



£ ,1, r „ Is ,1 

2 

so that we can take 

el 11 

(13) ev = -u 2 + n(a + -)u xt - efJ,[(b + -)uu xx + (c + -)u 2 ] 
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Step 3. We choose the coefficients j3, 7 and /j, such that the first equation of ((3]) is 
also satisfied up to 0(/i 2 ) terms. This is equivalent to checking that 

(14) [u + ev] t + [(1 + eu)u] x + e 2 [vu] x = 0(fi 2 ). 

First remark that one infers from (1131) that 



el 11 
ed t v = -uu t + n(a+ -)u xtt - efi[(b+ -)uu xx + (c+ -)u x ] t . 

= ~u(u x + e^uu x + ftauxxt) - (J,(a + ^)d xt (u x + e^(uu x )) 

+e^[{b + ^)uu xx + (c + \)u 2 x } x + 0(fi 2 ) 

1 2 3 2 

= -e-uu x - e -u u x - /Lt(a + -)u xxt 
5 5 

+e/i[(2a + b + -)uu xx + (-a + c+ l)u 2 x ] x + 0(p 2 ): 



similarly, one gets 

3 1 

e [vu] x = e 2 -u 2 u x - efi(a + - )[uu xx \ x + 0{^ 2 ), 

so that (|14p is equivalent to 

3 1 15 

u t + u x + e-uu x - [ie(a + -)u xxt = e\i [ - (a + b + -)uu xx - (-a + c + l)u 2 x ] + 0(/i 2 ). 

Equating the coefficients of this equation with those of (JT2J) shows that the first 
equation of is also satisfied at order 0(/j, 2 ) if the following relations hold: 

1,1 19 

a = , — , c= , 

6' 6' 48' 

and the conditions given in the statement of the proposition on a, (3, 7 and 5 follows 

from the expressions of a, b and c given after equation (|12p . □ 

As said in Remark [H none of the equations of the one parameter family con- 
sidered in Proposition [1] is completely integrable. This is the reason why we now 
want to derive a wider class of equations of the form ((7|) - and whose solution can 
still be used as the basis of an approximate solution of the Green-Naghdi equations 
(O (and thus of the water waves problem) . We can generalize Proposition Q] by 
replacing the vertically averaged velocity u given by ([2]) by the horizontal velocity 
u (9 & [0,1]) evaluated at the level line 9 of the fluid domain: 

so that 9 = and 9 = 1 correspond to the bottom and surface respectively The 
introduction of 9 allows us to derive an approximation consistent with (|3|) built on 
a two-parameter family of equations of the form . 

Proposition 2. Let p £ M, 9 £ [0, 1], and write X = \(9 2 — -|). Asssuming that 

, . 1 . 3 1 3, . 9 23 3, 

a = p + A, P = P h A, 7 = v A, = v A, 

6 2 6 2 2 24 2 

there exists D > such that: 

• For all s > and T > 0, 

• For all bounded family (u e ^ e ) (e ^ )eV € C([0, f ]; H S+D (R)) solving Q), 
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the family C e, ' 1 )(e.M)6'Pj with (ommiting the indexes £,n), 

(15) u = u e + y,\u B xx + 2/j,e\u e u e xx , 

(16) C : = U +7 W + ^u xt - £fJ,[- 



is consistent (of order s and on [0, ^\) with the Green-Naghdi equations f3J). 

Remark 2. i. The one parameter family of equations {?]) of Proposition]]^ corre- 
sponds to the particular case 9 2 = 1/3 (or A = 0J. 

ii. There exists only one set of coefficients such that 5 — 2"f and (3 — —27 (corre- 
sponding to p = — i and Q 2 = |, and thus a = — j, (3 = — < 0, 7 — 5 = 
The corresponding equation is therefore a Camassa-Holm equation in the sense of 
Definition 

3 15 5 

(17) u t + u x + -euu x - fi(-u xxx + —u xxt ) = —en(uu xxx + 2u x u xx ). 

iii. There exists only one set of coefficients such that 5 = 37 and (3 = — §7 (obtained 
with 6 2 = || and p = —-^fg, and thus a = — g|, (3 = 7 = gg, S = j^). The 
corresponding equation is therefore a Degasperis-Procesi equation in the sense of 
Definition [U 

3 .11 10 . 5 , 

u t + u x + -euu x - (i(—u xxx + —u xxt ) = —efi(uu xxx + iu x u xx ). 
2 54 27 36 

Proof. From the proof of Prop. 3.8 of [3J, one has the following expression (neglect- 
ing terms of order 0(/i 2 )): 

u = d x i> + M ((i + e()d x cd x ^ + (1 + 3 eC)2 dlip) , 

u 6 = d x if> + m((1 + eQd x Cd x il> + (1 + 2 £C)2 (1 - g 2 )d^), 

where -0 denotes the trace of the velocity potential at the surface. It follows from 
these formulas that 

= u e + t ^(9 2 l)u e xx + ve(6 2 - \)u e u e xx + 0(v 2 ), 
where we used £ = u e + 0(e) for the last equality. 

This formula, together with Proposition [TJ easily yields the result. □ 

2.2. Equations on the surface elevation. Proceeding exactly as in the proof of 
Proposition [Tl one can prove that the family of equations 

3 3 3 

(18) (t+(x+2£((x~g£ 2 ( 2 (x + Jj;£ 3 ( 3 (x+H(a(xxx+l3(xxt) = £H{j((xxx+6~(x(xx) 

for the evolution of the surface elevation can be used to construct an approximate 
solution consistent with the Green-Naghdi equations: 

Proposition 3. Let ggR and asssume that 

R 1 3 1 X 9 5 

a = o, p = q , 7 = — q , — — a . 

Then there exists D > such that: 
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• For all s > and T > 0, 

Tl. TTS + D, 



For all bounded family (C e ^) (e ^ )eV £ C([0, H s+D (R)) solving {7 
the family (C e,A \ u£,>1 )(e.^)e'P : with (omitting the indexes e, fx) 

"-■= c+ l(-i?-i< 3+ ^'- *k* + + 5<a) ■ 

is consistent (of order s and on [0, wi£/i t/ie Green-Naghdi equations J3J). 
Remark 3. Choosing q = 1/12, i/ie equation t!8\) reads 

(19) Ct + Cx + \<Cx ~ ^ 2 C 2 G + If 3 C 3 C. + -^(Cxxx - Cart) 

7 

While for any q G M, is an equation for the evolution of the free surface 

and all these equations have the same order of accuracy 0(e 4 ,fi 2 ), it is more 
advantageous to use U9\) since it presents better structural properties that we take 
advantage of in Q3.S\ The ratio 2 : 1 between the coefficients of CxCxx and CCxxx is 
crucial in our considerations. 

3. Mathematical analysis of the models and rigorous justification 

3.1. Large time well-posedness of the unidirectional equations ([7]) and 

(|18|). We prove here the well posedness of the general class of equations 

(20) u t +u x + \euu x + e 2 iu 2 u x + s 3 ku 3 u x + [i(au xxx + (3u xxt ) 

= e[i{juu xxx + Su x u xx ), 

with i, k 6 1; in particular, ([20]) coincides with (|7|) and (|18| if one takes i = k = 
and i = -f , « = 4 respectively. That is, we solve the initial value problem 



(21) 



u t +u x + \euu x + e 2 lu 2 u x + e 3 ku 3 u x + fJ,(au xxx + (iu xxt ) 



2 

= en(juu xxx + Su x u xx ), 

u U=a = u ° 



on a time scale 0(1/ e), and under the condition (3 < 0. In order to state the result, 
we need to define the spaces X s as 

Vs > 0, X s+1 = H S+1 (R) endowed with the norm \ f\ 2 x * + i = \ f\ 2 H s + fJ-\d x f\ 2 Hs , 
and we also recall that the set V is defined in (fTOl). 



Proposition 4. Assume that (3 < and let fi > 0, M > 0, s > | and u° £ 
H S+1 (M). Then, there exists T > and a unique family of solutions (ue^)( E .^)ev 
to fUl bounded in C([0, ^];X S+1 (1)) n ^([O, f]; X s 

Proof. For all v smooth enough, let us define the "linearized" operator C(v,d) as 
C(v, d) = (1 + n0d£)dt + d x + /.tad 3 + -evd x + e 2 iv 2 d x + e 3 nv 3 d x 

-efx-yvd 3 - e^8(^v x d 2 x + ^v xx d x ), 



In order to construct a solution to (|20[) by an iterative scheme, we are led to study 
the initial value problem 

f C(v,d)u = ef, 
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If v is smooth enough, it is completely standard to check that for all s > 0, / G 
L} oc (Wf; H S (R X )) and u° e H S (M), there exists a unique solution u G C(R+; (R)) 
to (|22j) (recall that /3 < 0). We take for granted the existence of a solution to (|2"2"1) 
and establish some precise energy estimates on the solution. In order to do so, let 
us define the "energy" norm 

Vs>0, E s (u) 2 = \u\\. - ii!3\d x u\ 2 HB . 

Differentiating ^e~ eXt E s (u) with respect to time, one gets, using the equation (|2"2"|) 
and integrating by parts, 

l -e M d t {e^ xt E a {u) 2 ) = - e -^-E s {u) 2 + e{K s f,K s u)-e{k s {Vd x u),k s u) 

+e^(A s (vd 3 u),A s u) - en-(A s (v x d x u),A s d x u), 

with V — \v + eiv 2 + e 2 nv 3 . 

Since for all constant coefficient skewsymmetric differential polynomial P (that is, 
P* = — P), and all h smooth enough, one has 

(A s (hPu), A s u) = ([A s ,h]Pu,A s u) - -([P,h]A s u,A s u), 

we deduce (applying this identity with P = d x and P = d x ), 

^e eXt d t {er eXt E s {u) 2 ) = -— E s {u) 2 - e([A s , V]d x u, A s u) + -((d x V)A s u,A s u) 
3 

-£/^y([A s , v]d 2 u - -v x A s d x u — v xx A s u,A s d x u) - efij([A s , v x ]d x u, A s u) 

-e^(A s (v x d x u), A s d x u) + e(A s f, A s u), 

Here we also used the identities 

[A s ,v]d 3 x u = d x ([A s ,v]d 2 x u) - [A s ,v x ]d 2 u 

and 

-(v xxx A s u,A s u) = -(v xx A s u,A s u x ). 

Since \u\h' < E s (u) and y/Ji\d x u\H' < ^_p E s {u)^ one gets directly by the Cauchy- 
Schwartz inequality, 

e EXt d t (e- EXt E s (u) 2 ) < eC(fi ,^,-y,6){A(u,v)E s (u) + B(v)E a (u) 2 ) 

P 

-eXE s (u) 2 + 2\f\ H sE s (u), 

with 

A(u,v) = \[A s ,V]d x u\ 2 + \[A s 1 v]d x (^JId x u)\ 2 + \[A s ,^JIv x }d x (^JId x u)\ 2 

+ y/Jj-\Vxd x u\ H -, 

B{v) = \d x V\ 

CO ~T | | OC 

+ \d x (y/JId x v)\ 

Recalling that for all s > 3/2, and all F, U smooth enough, one has 

|[A s ,F]C/| 2 <Cst \F\ H s\U\ H »-r, 
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it is easy to check that one gets A(u, v) < C(/to, M, Jj, I, K, E s (v))E s (u) and B(v) < 
C(no,M,±,i,n,E s {v)). Therefore, we obtain 

ef xt d t (e^ xt E s (u) 2 ) < (C(fx ,M, ± 7 , 6, t, «, £ s («)) - A)£» 2 + 2 £ ^(/)i? s ( t i). 

Taking A = At large enough (how large depending on C(/zo, M, 4,7, <5, t, k, £7 s (v))) 
to have the first term of the right hand side negative for all t £ [0, ~], one deduces 

V* G [0, j], d t (e- eXTt E s (u) 2 ) < 2ee- eXTt E s (f)E s (u). 
Integrating this differential inequality yields therefore 

Vf G [0, -], E s (u)(t) < e £XTt E°(u°) + 2e f e ^(t-*') E s {f{t'))dt' . 
e Jo 

Thanks to this energy estimate, one can conclude classically (see e.g. [1]) to the 
existence of 

T = T(fjto,M, \u°\ x s+i,-,-f,5,i,K,) > 0, 

and of a unique solution u £ C([0, ^]; X 8+1 (R d )) to (|21[) as a limit of the iterative 
scheme 

f C(u n ,d)u n+1 = 0, 

u = u , and Vn G N, 4 «+i _ „o 

til Ct . 

^ |t = 

Since u solves (f2"0|. we have C(u,du)u = and therefore 

(A s_1 (l + nPd%)dtu, A'^dtu) = -e{h s - l M{u,d)u, A" -1 ^), 
with .M(u, 9) = £(w, 9) — (1 + iif3d 2 )d t - Proceeding as above, one gets 

1 

and it follows that the family of solution is also bounded in C 1 ([0, — ]; X s ). □ 



E'-^Btu) < C(jm, M, \u°\ x s+i, -, 7 , 6, i, k, £ s (u)), 



3.2. Rigorous justification of the unidirectional approximations ([7J. In 

Proposition^ we constructed a family (u e,A1 , £ £ >^) consistent with the Grccn-Naghdi 
equations in the sense of Definition |4] A consequence of the following theorem is 
a stronger result: this family provides a good approximation of the exact solutions 
(w £ ' M , C e,AI ) of the Green-Naghdi equations with same initial data in the sense that 
(M^C ) = (« e ' At ,C E '' i ) + 0(^ 2 t) for times 0(l/e). 

Theorem 1. Let /io > 0, M > 0, T > and "P be as defined in U0\) . Let also 

p £ R ; G [0, 1], and a, (3, 7 and <5 6e as in Proposition^ 

If(3<0 then there exists D > and T > sucft fftaf /or a// it G H S+D+1 (R): 

• There is a unique family (u £ ^ 7 C'^Cs^^eV G C([0, H s+D (R) 2 ) given by 
the resolution of with initial condition u° and formulas jl5\) - [Tfy) : 

• There is a unique family (u Slft ,C ,fl )(e,pi,)eV <S C([0, j}; H s+D (M.) 2 ) solving 
the Green-Naghdi equations f3|) with initial condition (u e,A1 , C e ''')|t=o ■ 

Moreover, one has for all (e, /1) £ V , 

Vt £ [0, |], |n e ^ - u e '"U- ([0jt]xH ) + - ^'"U-ao.tixR) < Cst M 2 i. 
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Remark 4. It is known (see that the Green-Naghdi equations give, under the 
scaling a correct approximation of the exact solutions of the full water waves 
equations (with a precision 0([i 2 t) and over a time scale 0(\/e)). It follows that 
that the unidirectional approximation discussed above approximates the solution of 
the water waves equations with the same accuracy. 

Remark 5. We used the unidirectional equations derived on the velocity as the basis 
for the approximation justified in the theorem. One could of course use instead the 
unidirectional approximation 118\) derived on the surface elevation. 

Proof. The first point of the theorem is a direct consequence of Proposition 2] 
Thanks to Proposition [2j we now that (u e,,1 ,( £ ' fJ ')e,n is consistent with the Green- 
Naghdi equations ([3]), so that the second point of the theorem and the error estimate 
follow at once from the wcll-posedness and stability of the Green-Naghdi equations 
(see Th. 3 of [4]- note that instead of using this general result which holds for two 
dimensional surfaces and nonflat bottoms, one could easily adapt the simpler and 
more precise results of [25] to the present scaling). □ 

3.3. Wave breaking. For the Camassa-Holm family of equations ((7]) for the ve- 
locity it is known (see [8]) that singularities can develop in finite time for a smooth 
initial data only in the form of wave breaking. We will show now that this form of 
blow-up is also a feature of the equation (jTSJ) for the free surface. More precisely, if 
a smooth initial profile fails to produce a wave that exists for all subsequent times, 
then we encounter wave breaking in the form of surging (and not plunging, as would 
be the case if ([7]) were the equation for the evolution of the free surface). 

Our first result describes the precise blow-up pattern for the equation (|19|) for 
the free surface. 

Proposition 5. Let Co G H 3 (M.). If the maximal existence time T > of the 
solution of U9\) with initial profile C(0, •) = Co * s finite, T < oo, then the solution 
(E C* 1 ([0,T);i7 2 (R))nC([0,T);i7(]R)) is such that 

(23) sup {|C(*,x)|} <oo 

te[0,T), x£R 

and 

(24) sup {&:(*, a;)} T oo as 1 1 T. 

Proof. In view of Proposition |4] given Co S H 3 (R) , the maximal existence time of 
the solution to (fl9|) with initial data C(0) = Co is finite if and only if |C(^)|ff 3 (R) 
blows-up in finite time. Thus if (|24|) holds for some finite T > 0, then the maximal 
existence time is finite. To complete the proof it suffices to show that 

(i) the solution ((t) given by Proposition [4] remains uniformly bounded as long 
as it is defined; 

and 

(ii) if we can find some M = M(Co) > such that 

(25) Cx(t,x)<M, lei, 

as long as the solution is defined, then |C(0l-ff 3 (K) stays bounded on bounded time- 
intervals. 
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Item (i) follows at once from the imbedding L°°(R) C H 1 (R) since multiplying 
l9|) by £ and integrating on M yields 



15 s J CCxxCxxx d 8 / C CxxCxxx dx 



(26) 8 t ( /[C 2 + ^M Ul\dx) =0. 

y JR lz JR ' 

To prove item (ii), notice that multiplication of (|19|) by Caraa; yields 

(27) dt( [ [CL + ^M / cLJ^) 

+ T6 £3 I ^ dx + Y £3 I dx+^fie ( CxCxx dx. 

after performing several integrations by parts. If (|25[) holds, let in accordance with 
(f2"6"| the constant Af > be such that 

|C(t,a?)| < Af 0j xeR, 

for as long as the solution exists. Using the Cauchy-Schwartz inequality as well as 
the fact that fi < 1, we infer from (|26p and ([27)) that 



„ , s /90e 45e 2 , 27e 3 , 45e 3 , 

S t E(t) < ( — Mo + — M 2 + — M 3 + — M 3 + 21 eM ) £(*), 



where 



= ^IC 2 + ^ K 2 + CL + ^ MCLJ da;- 
An application of Gronwall's inequality enables us to conclude. □ 



Our next aim is to show that there arc solutions to (|19[) that blow-up in finite 
time as surging breakers, that is, following the pattern given in Proposition [5] We 
will prove this by analyzing the equation that describes the evolution of 

(28) M{t)= sup {CzQt,x)}. 

xeR 

For the degree of smoothness of the solution ((t) given by Proposition [U we know 
that M(t) is locally Lipschitz and 

(29) ± M (t) = ( tx (t,at)) fora.e. t, 

where £(i) is any point where M(t) = Cx(t,£,(t)) cf. [5]. For further use, let us also 
note that 

(so) (i-±-nd 2 x )- 1 f = p*f, /a 2 (i), 



where 
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^3 



3 -2j-\x\ 



with 

(3D 11^- = ^, 11*11* = 1, 

and 

(32) \\P x \\l~ = ~, \\Px\\li = 2 J^-, \\Px\\l^V2(-Y <4^i. 
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Applying (1 — jk \l d 2 )" 1 to (fTO|) . we obtain the equation 

Ct + P X *C + \eP x *( 2 - \s 2 P x *C 3 + |r£ 3 ^ * C 4 
4 8 64 

Differentiating this equation with respect to the spatial variable, we obtain 
C tx +dtP*(+ledlP*e- 1 -e 2 d 2 P*e + ^d 2 x P*C* 

+ ^^P*C = ~I^ d 2 x P *C-lf,eP x * (CCxxx). 
Since CCxxx = <9 2 (CG0 - 3( x ( xx and 

(33) d 2 x P*f = P x *t x = —P*f-—f, fEL 2 (R), 

li fx 

we deduce that 

(34) ( tx + 2P X * ( x -\e 2 P x * (C 2 C.) + % Px * (C 3 C.) 

8 fo 

= U + 7 {p*£ + 7 -^£ + y CC*, - ^p, * (CCx). 

We can now prove the following blow-up result. 
Proposition 6. If the initial wave profile Co G H 3 (R) satisfies 

SU P {Co(x)} 2 > 2 l Co ^ + - \eCl' 2 ^"3/4+1^^-3/4 
, 7 r ,-1/2 , 16 r l/2 -3/4.-1 

+ 2°oM ~3~ ^ ' 

w/iere 

C =^[Co + (Co) 2 ]^>0, 
i/ien wave breaking occurs for the solution of U9\) m finite time T = O(^). 
Proof. Notice that 

sup{|C 2 (*)|}< J Ue + C)dx = % 



Therefore, using Young's inequality and the estimates (|31[) - (|32j) . we obtain that 

P**&\\l- < \\Px\\lA\U\l* <4 M - 3/4 C 1/2 , 
P*C 2 ||l~ < ||P|U-||^|Ui < \\P\\l-\\Q\1* < 2^ 1/2 c , 

p x *(CC*)\\l- < \\Px\\l4C(x\\li < IIPxIU'llCIU-HCrlU' <4/i- 3 / 4 c , 
Px*(C 2 <x)h°° < \\p x \\lA\C 2 CAl- < \\p x \\lA\C\\1~\\Kx\\l* <2v- 3/4 c* /2 , 
Px*(C 3 Cx)lk~ < \\Px\\l4C 3 CAl^ < ||Px||L 2 ||Clli=c.||||CllL 2 <2m" 3/4 C 2 - 



Since (|34|) is at any fixed time an equality in the space of continuous functions, we 
can evaluate both sides at some fixed time t at a point £ M where M(t) = 
Cx (*,£(*)), with M(*) defined in ([28]). Since ( xx (t,£(t)) = 0, from ([29]), and the 
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previous estimates we derive the following differential inequalities for the locally 
Lipschitz function M(t): 

(35) \ e M 2 (t) + (14 Co e + \ e 2 C 3/2 + ^ C 2 + 8C 1/2 ) M " 3 / 4 + \ eC^' 2 
> M'(t) for a.e. t, 



and 



(36) M'(t) > L e M 2 (t) - (14 C e + | e 2 C 3/2 + | e 3 C 2 + 8C 1/2 ) ^ 

for a.e. t. 

Notice that ^ Co G H 3 (R) ensures M(0) > 0. At i = the right-hand side of 

([3"6"| is by our assumption on the initial wave profile larger than -eM 2 (0). Wc 

infer that up to the maximal existence time T > of the solution of (fTl?)) the 
function Af(t) must be increasing and, moreover, 

M'(t)>-eM 2 (t) for a.e. i. 

Dividing by M 2 (t) > M 2 (0) > and integrating, we get 

111 r 

>-ei, <e[0,T). 



M(0) M(t) ~ 4 
Therefore lim nT M (t) = oo and T < 



eM(0)' 

On the other hand, a similar argumentation applied to (|35j) yields 
M'(i) < AeM 2 {t) for a.e. i, 

so that 

1 1 

— rr > — tt - 4st 
M(t) ~ M(0) 

as long as the solution of (fT9|) is defined. Since lim t |t M (t) = oo wc deduce from 

the previous inequality that T > — . Thus the finite maximal existence time 

4eM(0) 

T > is of order 0(±). □ 

3.4. Numerical computations. In this section, we use numerical computations 
to check that the surface equation (|19[) and the Camassa-Holm equation (fT7|) lead 
respectively to surging and plunging breakers as predicted theoretically by Propo- 
sition [6] (and [8] for (fTTjl ). We use the same kind of numerical scheme for both 
equations; in fact, our scheme works for any equation of the class (|20| with f3 < 0. 
In order to reduce the size of the computational domain, we solve (|20|) in a frame 
moving at speed 1, in which (|20p is replaced by 

Ut + f^uuj; + £ 2 iu 2 U£ + e 3 ku 3 U£ + //((a — /3)u£££ + (3u^t) 
where £ stands for a; — t. 

The numerical scheme used here is a simple finite difference leapfrog/ Crank-Nicolson 
scheme whose scmi-discrctized version reads 

,,n+i _ ,,n— l 
Vn > 1, (1 + m/39 2 ) ^ = 
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with 

3 

F[u] = —-euu£ — e 2 iu 2 u£ — e 3 ku 3 u£ — fj,(a — /3)w^ + e/-i( , -fuu^ + Su^u^), 

and where St is the time step and u n ~ U\ t=nS (to start the induction, that is, for 
ri = 0, the centered discrete time derivative must be replaced by an upwind one). 
Numerical computations are performed for (|1T[1 and (|19|) with the same initial value 

U| t=0 = cxp(-100a; 2 ), 

and with fi = 0.2, e = y/Ji. Figure [T] shows the formation of a plunging breaker for 
the solution of (| 1 T[) : the little mark on the curves materializes the point of minimal 
slope. For the same initial data, Figure [2] shows the formation of a surging breaker 
for the solution of (fl9|) ; the little mark on the curves materializes here the point of 
maximal slope. 




FIGURE 1. A plunging breaker for the Camassa-Holm equation (fT"7|) 
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